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THE BEHAVIOUR OF THE COMPLETE EIGENSTRUCTURE OF A 
POLYNOMIAL MATRIX UNDER A GENERIC RATIONAL 
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Abstract. Given a polynomial matrix P(x) of grade g and a rational function x(y) = n{y)/d{y), 
where n{y) and d(y) are coprime nonzero scalar polynomials, the polynomial matrix Q{y) := 
[d{y)]s P{x{y)) is defined. The complete eigenstructures of P{x) and Q{y) are related, including 
characteristic values, elementary divisors and minimal indices. A Theorem on the matter, valid in 
the most general hypotheses, is stated and proved. 
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1. Introduction. A polynomial matrix is a matrix whose entries belong to some 
polynomial ring i? H] . In this paper we will always assume that i? is a principal ideal 
domain. This condition is equivalent to i? = F[x], the ring of univariate polynomials 
in X with coefficients lying in some field F. 

An important property of a polynomial matrix with entries in ¥[x] is its complete 
eigenstructure, whose definition is given in Subsection 12.21 The name eigenstructure 
comes from the special case where F = C; in this context, polynomial matrices are 
usually seen instead as matrix polynomials, that is polynomials whose coefficients 
are matrices [?]■ Any matrix polynomial is associated with a polynomial eigenvalue 
problem (PEP); the complete eigenstructure is strictly related with the properties 
of the associated PEP. More precisely, it gives the complete information about the 
eigenvalues, eigenvectors and Jordan chains of the matrix polynomial, and also about 
the Kronecker form of any strong linearization of the matrix polynomial [T]. Poly- 
nomial eigenvalue problems arise in many applications, from mathematics, science 
and engineering; both their algebraic properties and the numerical methods for their 
approximate solutions are widely studied. See, e.g., [51 [T^ [TB] . 

The aim of this paper is to investigate the link between the complete eigenstruc- 
tures of two polynomial matrices P{x) and Q{y) related one to another by a rational 
transformation x{y) of the variable. In order to better explain the question we are 
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interested in, let us consider the following example, where R — 
we have to deal with the polynomial matrix 



Suppose that 
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if we choose grade(P(a;))= 2 (the grade of P{x) is an arbitrary integer g such that 
g > degP{x); more details are given in Section [2]), then the complete eigenstructure 
of P{x) is the following: 

• the elementary divisors of P{x) are {x — 20), x, {x — 20), x'^; 

• there are no right minimal indices; 

• the left minimal indices of P(x) are 0, 1. 

The rational change of variable x{y) = ^^^iZy^ induces an application $2, as defined 

in (EU, such that <f2{P{x)) - (y^ - yfPC-^f^) =: Q{y), with grade(Q(y)) = 4 
(see Section [3]) and 



Q{v) 



(25-162;2)(2y-5)2 
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By studying the complete eigenstructure of Q{y) we find out that 



• the elementary divisors of Q{y) are (y 

• there are no right minimal indices; 

• the left minimal indices of Q{y) are 0, 2 



i)2 



Notice that x(|) = 20, x{±j) = 0, and that y = | is a root of multiplicity 2 of 
the equation x{y) = 20 while y = ±| are roots of multiplicity 1 of the equation 
x{y) = 0. We can therefore conjecture that if {x — xqY is an elementary divisor of 
P{x) and yo is a root of multiplicity m of the equation x{y) = xq then (y — yo)™'^ is 
an elementary divisor of Q{y). Moreover, we see that apparently the minimal indices 
have been multiplied by a factor 2; notice that 2 is the degree of the considered 
rational transformation (that is the maximum of the degrees of the numerator and 
the denominator). 



The main result of the present paper is the proof that the conjectures above, which 
will be stated more precisely in SectionlH are true for every rational transformation of 
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the variable x{y) and every polynomial matrix P{x). Moreover, analogous properties 
hold for infinite elementary divisors and right minimal indices. 

The motivation for this work comes from the will to generalise the partial results 
derived in [5] , where we considered the particular case of a square and regular poly- 
nomial matrix with entries in C[a;] and without infinite elementary divisors, and the 

2_|_-| 

Dickson change of variable x{y) = ^ J" . Moreover, we wish to extend the results by 
D. S. Mackey and N. Mackey [H], who described the special case of rational transfor- 
mations of degree 1, also known as Mobius transformations. The present contribution 
is offered as both a synthesis and an extension of the previous works cited above. 

The results provided in this paper can be used to design numerical methods for 
the approximate solution of PEPs. An example in this regard, restricted to the case 
of the Dickson transformation, is given in [5] for the solution of the palindromic PEP. 

The structure of this paper is the following: in Section[2]we expose the theoretical 
background we are going to work within, and we give some basic definitions that we 
will use later on. In Section |3] we formally define the application between polynomial 
matrices induced by a rational change of variable and we present some intermediate 
results. Our main result is Theorem 14.11 which is stated and commented in Section 
m Sections [5] and m are devoted to the proof of our result. For the sake of simplicity, 
in Sections [5] and [5] we assume that the underlying field is algebraically closed: in 
Section [7] we show how the result still holds for an arbitrary field. Finally, in Section 
[51 root polynomials are introduced in order to prove a technical Lemma. 

The first part of Theorem 14.11 was stated and proved, but only for a very special 
case, in [5]. Besides the generalisation to a generic rational transformation and a 
generic polynomial matrix, this paper also contains the analysis of what happens to 
minimal indices and infinite elementary divisors. 

2. Preliminary definitions. In this Section we describe our notation and recall 
some basic definitions. 

2.1. Basic facts on polynomials. Let Z he a ring and let Z[x] be the ring of 
the univariate polynomials in the variable x with coefficients in Z. We denote the 
degree of z S Z[x] by the letter k, and sometimes write k — degz. 

On the other hand, the grade |10j of a polynomial z e Z[x] is any integer g = 
grade(0) satisfying g > k. The choice of the grade of a polynomial is arbitrary: 
nevertheless, some algebraic properties of polynomial matrices depend on the grade. 

Remark 2.1. In some sense, the degree of a polynomial is an intrinsic property 
while the grade depends on its representation. In fact, informally speaking, the grade 
depends on how many zero coefficients one wishes to add in front of the polynomial. 



4 



V. Noferini 



Let now g be the grade of z = X]i=o '^i^^ ^ Z[x\. The reversal of z with respect 
to its grade O [10] is 



1=0 

The subscript 5 will sometimes be omitted when the reversal is taken with respect to 
the degree of the polynomial, that is RevfeZ =: Revz. 

Let now F be an arbitrary algebraically closed field. 

Remark 2.2. Although the hypothesis that F is algebraically closed is useful to 
state in a simpler way our results, it is not strictly necessary. See Section^ 

A well-known result that is crucial to us is that F[a;] is guaranteed to be an 
Euclidean domain. Given zi, Z2 G ¥[x], not both zero, we denote by GCD(zi, Z2) their 
greatest common divisor; we additionally require that GCD(zi, Z2) is always monic so 
that it is uniquely defined. We say that zi and Z2 are coprime if GCD(zi, Z2) = lF[a;]. 

Notice that a polynomial z G F[x] can be thought of as a function z{x) : F — > F. 
Thus, applying (|2.ip . in this case the formula Revgz(x) — x^z{x^^) holds. 

Let now Z™^p be the set of ni x p matrices with entries in Z; the case p = 1 
corresponds to the set of vector with m elements in Z, denoted by Z™. We are mainly 
interested in analysing (F[a;])'"^^, the set of m x p polynomial matrices with entries in 
F[a;]. Mm(F[a;]) := (F[x])™'*™ is the ring of square polynomial matrices of dimension 
m. A square polynomial matrix A e Afm(F[a;]) is said to be regular if det A 7^ Ovix] 
and singular otherwise. If A is regular and det A CzV then A is called unimodular. 

Remark 2.3. Notice that {¥[x])'"^^p = {¥^'^^p)[x]; or in other words, a poly- 
nomial matrix, defined as a matrix whose entries are polynomials, is also a matrix 
polynomial, defined as a polynomial whose coefficients are matrices. 

The notions of grade and degree can be extended in a straightforward way to 
polynomial matrices, as follows: the grade (resp., the degree) of A e {Z[x])"^^p is 
defined as max.ijgrade(ylij) (resp., as maxij degAij). Analogously, the reversal of a 
polynomial matrix is defined just as in (|2.ip . after replacing ai ^ Z with Bi € Z"^^p. 

2.2. Characteristic values, elementary divisors, and minimal indices. 

Let A e {F[x])"^^P , and let 1/ =: mm{m,p). Suppose that there exist Di, . . . ,Di, G 
¥[x] such that Aij — DiSij, where Sij is the Kronecker's delta. Then we write A = 
diag(Di, . . . , D,^), and we say that A is diagonal. Notice that we use the notation 
indifferently for both square and rectangular polynomial matrices. 

The following Theorem, which in its most general is due to Frobenius [3], is in 
point of fact valid for any matrix with entries in a principal ideal domain [4l [7] . 
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Theoreim 2.4. Let P{x) e (F [a;] Then there exist two unimodular A{x) S 
Mjn{^[x\) and B{x) e Mp{¥[x]) such that 

S{x) = A{x)P{x)B{x) = diag(di(x), . . . , d^{x)), 

where di{x) G F[a;] is monic ^i < v :^ niin(m,p) and di{x)\di+i{x) Vi < — 1. 

Such an Six) e (F[a;])™^P is called the Smith form [a[14j of and the 

are called its invariant polynomials [4j . The Smith form, and thus the invariant 
polynomials, are uniquely determined by P{x). Notice that a square polynomial 
matrix P{x) is singular if and only if at least one of its invariant polynomials is zero. 

Using the fact that F is algebraically closed, let us consider a factorization of the 
invariant polynomials over F[x]: di{x) = Yiji^ ~ a^j)*^^''"'- The factors {x — a;_,)'^^'<'' 
are called the elementary divisors of P{x) [H [7| corresponding to the characteristic 
value Xj [4]. Notice that, from Theorem 12. 4| ii < 12 ^ %,(«2)- 

Remark 2.5. When F = C the characteristic values of the polynomial matrix 
P{x) are often called the eigenvalues of the matrix polynomial P{x). Given an eigen- 
value xq, there is a Jordan chain of length £ at xq if and only if {x — xqY is an 
elementary divisor. The number of Jordan chains at xq is equal to the number of 
invariant polynomials that have Xq as a root J^. 

Let us now denote by F(x) the field of fractions of the ring F[a;]. Let V be a vector 
subspace of (F(a:))^', with dimV = s. Let {vi} be a polynomial basis for V with the 
property deg vi < ■ ■ ■ < deg Vs . Often we will arrange a polynomial basis in the 
matrix form V{x) = [vi{x), . . . ,Vs{x)] S (F[x])^^*'. Clearly, polynomial bases always 
exist, because one may start from any basis with elements in the (vectorial) field of 
fractions, and then build a polynomial basis just by multiplying by the least common 
denominator. Let := deg Vi be the degrees of the vectors of such a polynomial basis; 
the order of V{x) is defined [2] as '^i- ^ polynomial basis is called minimal [2] if 

its order is minimal amongst all the polynomial bases for V, and the are called its 
minimal indices [2]. It is possible to prove [21 14] that, although there is not a unique 
minimal basis, the minimal indices are uniquely determined by V. 

The right minimal indices [Tj of a polynomial matrix P{x) € {¥[x])"^^^ are defined 
as the minimal indices of kerP(a;). Analogously, the left minimal indices [I] of P{x) 
are the minimal indices of kerP(a;)^. 

Given the grade g of P{x), we say that 00 is a characteristic value of P{x) if Op 
is a characteristic value of ReVgP{x). The elementary divisors corresponding to 00 
are defined [9] as the elementary divisors of ReVgP{x) corresponding to Of; if x^ is 
an elementary divisor of ReVgP{x) we formally write that {x — ooY is an infinite ele- 
mentary divisor of P{x). Notice that the infinite elementary divisors of a polynomial 
matrix clearly depend on the arbitrary choice of its grade. 
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We complete this section with the following definition pj: the complete eigen- 
structure of P{x) is the set of both finite and infinite elementary divisors of P{x) and 
of its left and right minimal indices. 

3. Rational transformations of polynomial matrices. Let n{y)^ d{y) G F[y] 
be two nonzero, coprime polynomials. Let us define N := degn(y), D := degd(y), 
and G :— ma.x{N, D). We will always suppose G > 1, that is n{y) and d{y) are not 
both elements of F. We denote the coefficients of n{y) and d{y) as Ui G ¥, i = 0, N 
and e F, j = 0, . . . , D, that is n{y) ^ J2fLo ^-iv) = J2Za '^^v' ■ 



-1 



Let us introduce the notation F* := F U {oo}, having formally defined oo := Op . 
We consider the generic rational function from F* to F*: 

The function (13. ip induces an application ^g^n{y),d(y) ■ (F[a;])™^P (F[y])™^P defined 
as 

%MyUiy) (Pi^)) = Q{y) [d{y)rPix{y)) (3.2) 

Here g is the grade of P{x) G {¥[x])™^p, so for any choice oi g a different application 
is defined. We will usually omit the functional dependence of $ on n{y) and d{y) 
unless the context allows any possible ambiguity; also, if the grade is chosen to be 
g — k we will sometimes omit the subscript g, that is <^{P{x)) := ^k.n(y).d(y)iP{^))- 

Since a polynomial matrix is also a matrix polynomial, we can write P{x) = 
J2i=o Pi^^ ^'-'^ some Pi G F'"^*', i = 0, . . . , g. Notice that following the same point of 
view we can also write Q{y) — X)f=o 

Lemma 3.1. degQ{y) — deg$g(P(a;)) is less than or equal to q :— gD + 
niaxi:P;^o(*^ — iD)- If N D the strict equality degQijj) = q always holds. More- 
over, q < gG. 

Proof. Writing Q{y) as above, we can see it as a sum of the fc + 1 polynomial 
matrices Qi{y) = Pi[n{y)Y[d{y)]^~^ , < i < k, with either Qi{y) — Pi — oi 
deg Qi{y) = gD + i{N — D). Since the degree of the sum of two polynomials cannot 
exceed the greatest of the degrees of the considered polynomials, degQ(y) cannot be 
greater than q. Notice that if = G then gG > q = kG+{g — k)D and the maximum 
is realised hy i = k, while otherwise the maximum is realised by the smallest index 
j such that Pj ^ 0, and q = {g - j)G + jN. This means that if < G and Po = 
then q < gG, while q = gG if N < G but Pq ^ 0. 

Notice finally that, if ii ^ 12, then Qi-^{y) and Qi2(jj) have the same degree 
if and only if D = N . Since degQi,(y) ^ degQ^Jy) ^ deg{Qt,{y) + Qi^iy)) = 
ma.x{degQi^{y) , degQi.^{y)) , D ^ N is a, sufficient condition for deg Q{y) = q. D 
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Lemma [3.11 shows that degQ(y) < 9 < gG. The next Proposition describes the 
conditions under which the equahty degQ{y) — gG holds. 

Proposition 3.2. Let Q{y) ~ $(P(a;)). It always holds degQ{y) < gG, and 
degQijj) < gG if and only if one of the following is true: 

1. N > D and g > k; 

2. N < D , and there exist a natural number a > 1 and a polynomial matrix 
P{x) £ (F[a;])'"^P such that P{x) (x - x)''P{x), where x := ncd^^ if 
N = D^Gandx:=QrifN<D = G. 

Proof. Lemma [3.11 guarantees degQ{y) < gG. To complete the proof, there are 
three possible cases to be analysed. 

• If G = iV > 13, we know from Lemma lOI that deg Q{y) — q, and in this case 
q = gD + kN - kD. Therefore, degQ(y) = gG g ^ k. 

• If iV = D = G and, we get q = gG. Let Q{y) = Yl=o ^if- then, deg Q(y) < 
gG 4^ 9gG = 0(F[a;])mxp. On the other hand Qgc is the coefBcient of y^'^ in 

Qiy) = ELo P^Hy)my)V-\ so e,G = 4 ELo ^^^g-^/ = d%p{nGd-^^). 

Therefore, Qgc is zero if and only if every entry of P{nQd(^) is equal to Q-w[x]^ 
or in other words if and only if P{x) = {x — nQd^)"^ P{x) for some a > 1 and 
some suitable polynomial matrix P{x). 

• li N < D = G^ recalling the proof of Lemma [5TT] we conclude that deg Q{y) < 
gG if and only if Pq = 0, which is equivalent to P{x) = x°-P{x) for a suitable 
value of a > 1 and some polynomial matrix P{x). 

□ 

The grade of Q{y) is of course arbitrary, even though it must be greater than 
or equal to its degree. Since degQ(y) < q < gG, we shall define that the grade of 
Q{y) is gG. This choice has an influence on the infinite elementary divisors of Q{y), 
as they are equal to the elementary divisors corresponding to zero of the reversal of 
Q{y) taken with respect to its grade, that is Rev(gQ)(5(y). 

If one is interested in picking a different choice for the grade of Q{y), the following 
Proposition explains how the infinite elementary divisors change. 

Proposition 3.3. Let P{x) e {¥[x]Y'''^p, with k = degP(x). Then the finite ele- 
mentary divisors and the minimal indices of P{x) do not depend on its grade, while the 
infinite elementary divisors do. Namely, let v = min(m,p); x^^^di{x), . . . ,x^^^di^{x) 
are the invariant polynomials o/ReVgP(a;) if and only if di(x), . . . ,du{x) are the in- 
variant polynomials ofKevkPix), for any choice of g> k. 

Proof. Neither Theorem 12.41 nor the properties of kerP(x) and kerP-^(a;) de- 
pend on the grade, so minimal indices and finite elementary divisors cannot be af- 
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fected by different choices. Let S{x) — A{x)KevkP{x)B{x) be the Smith form of 
RevkP{x). We have RevgP(x) = a;3~'^RevfcP(x), which imphes that x^~'^S{x) = 
A{x)RevgP{x)B{x). Clearly di{x)\dj{x) xS^''di{x)\x^~''dj{x), and therefore we 
conclude that x^~'^S{x) is the Smith form of RcvgP{x). □ 

Let a, /3, 7, (5 e F. If G = 1, '^g^ay+p,'yy+s is clearly invertible and its inverse, with 
a little abuse of notation, is ^g^p-Sx^^yx-aiQiy)) = [^x — a]^Q( ;^~^^ ) . The most 
general case is analysed below. 

Proposition 3.4. Let us denote by F[x]g the set of the univariate polynomials 
in X whose degree is less than or equal to g. Given g,n{y),d{y), the application 
^g,n{y),d{y) '■ (■'I^'Wg)™^^ ~^ (^[y](gG))"^^^ always an injective function, but it is not 
surjective unless G — 1. 

Proof. Notice that $g can be thought as acting componentwise, sending P{x)ij 
to Q{y)ij = ^g{P{x)ij). Thus, it will be sufficient to show that, in the scalar case 
$g : ¥[x]g — > F[?;](gG), $g is surjective if and only if G = 1. This is true be- 
cause any polynomial that does not belong to the set Ry := {a{y) G ¥[y] : a{y) = 
Sf=o cannot belong to the image of $g, and Ry = F[y](g(3) if and 

only if G = 1. In fact, if we require that a generic r £ ¥[y]{gG) belongs to Ry, we find 
out that the g + 1 coefficients flj must satisfy gG + 1 linear constraints. 

To prove injectivity: ^g{Pi{x)) = ^g{P2{x)) ^ Pi{x{y)) = P2{x{y)) ^ Pi{x) = 
P2{x). □ 

Proposition 13.41 tells us that, unless G = 1 (the Mobius case), not every Q{y) is 
such that Q{y) — <f>(P(a;)) for some P{x). 

A couple of additional definitions will turn out to be useful in the following. Let 
xq G F*: we define T^q as the counterimage of a;o under the rational function x{y). 

Moreover let a, /3 G F be such that ^ = a;o and a and (3 are not both zero. For 
instance, we can pick (a,/?) = (a::o,lF) \i xq ^ oo and (a,/?) = (1f,0f) otherwise. 
Consider the polynomial equation 

ad[y) = Pn{y). (3.3) 

Let S be the degree of the polynomial adly) — j3n{y). Equation p.3p cannot have 
more than S finite roots. If S* < G then we formally say oo G Tx^. 

Remark 3.5. Notice that there are three cases that lead to S < G: 

1. N = D = G and xq — ncdQ^ , so that p.3p becomes dcniy) — nGd{y): in 
this case, S is the maximum value of i such that Ui ^ Xodi; 

2. N < D = G and xq = Of, so that p.3p becomes n{y) = Of and S — N ; 

3. D < N — G and xq = oo, so that p.3p is d{y) — Of and S — D. 
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We now define the multiplicity toq of any finite yo € T^g as the multiplicity of yo 
as a solution of the polynomial equation p.3|) . If oo € T^g , its multiplicity is defined 
to be equal to G — S. Therefore, the sum of the multiplicities of all the (both finite 
and infinite) elements of T^^ is always equal to G, while the sum of the multiplicities 
of all the finite elements of T^g is S. 

The finite elements of Txg are characterised by the following Proposition. 

Proposition 3.6. Let yo € F and xo e F*. Then yo e T^g if and only if 
yo is a solution of (j3.3|) for a,fi : xq — ^. Moreover, aid{yo) = /3in(yo) (ind 
a2d{yo) = /32n{yo) if and only if ^ ^ ^. 

Proof. The definition of T^g implies the first part of the Proposition. The second 
part comes from the fact that x{y) is a function. □ 

Proposition 13.61 albeit rather obvious, has the following important implication: 

Corollary 3.7. xo ^ xi -i^ Txg n T^^ = 0. Equivalently, ai/32 ^ Q.2P1 if and 
only if [I3in{y) — aid{y)\ and [/32".(y) — a2d{y)] G F[j/] are coprime. 

In particular, for any finite xq (z ¥, ^{x — xq) and d{y) are coprime. 

In order to clarify the latter definitions, let us consider an example. Let F = C 
and take n{y) — y"^ + y^ — y'^ — y + 1, d{y) = y^. Ti is the set of the solutions of the 
equation n{y) = d(y), so in this case Ti — { — 1,1, 00}. Moreover, the multiplicity of 
— 1 and 1 are, respectively, 1 and 2; since S — i and G = 4, the multiplicity of 00 is 
by definition G — S ~ 1. Within the same example. Too = {0}; has multiplicity 4 
because it is a root of order 4 of the equation d{y) — 0. 

4. Main result. We are now able to state our main Theorem. 

Theoreim 4.1. Given m^p g Nq andn{y),d{y) € ¥[y], let Xq e ¥* be a character- 
istic value of P{x) S (F[a;])"'^P, and let {x — XqY^, . . . , [x — XqY^ he the corresponding 
elementary divisors. Let g be the grade of P{x), define G — max(degn(y),degc?(?/)) 
and let gG he the grade of Q{y) = := [d(y)]9P(^) G (F[y])"><J'. Then for 

any yo G T^g : 

• yo is a characteristic value of Q{y); 

• (y — yo)™"^S • ■ • , (y — yo)"*"^^ o-re elementary divisors corresponding to yo 
for Q{y), where uiq is the multiplicity ofyo- 

Conversely, if Q{y) = ^g{P{x)) for some P{x), and if yo G F* is a characteristic 
value of Q{y) with corresponding elementary divisors (y — yo)"^ , • . • , (y ~ 2/0)"^ •' 

• xo — is a characteristic value of P{x); 

• mo\Ki yi < j, where mo is the multiplicity ofyo as an element ofTx„, and 
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{x — xo)"^o ''ij ...^ {x — Xq)"^" are elementary divisors corresponding to 
xq for P{x). 

In addition, the following properties hold: 

• the right minimal indices of P{x) are /3i, /3s if and only if the right 
minimal indices of Q{y) are G/3i, . . . , G/3s; 

• the left minimal indices of P{x) are 71, . . . , 7,- i/ and only if the left minimal 
indices of Q{y) are Gji, Gjr- 

For any choice of the apphcation $g, Theorem 14.11 gives a thorough description of 
the complete eigenstructure of $g(P(a;)) with respect to the complete eigenstructure 
of P{x). Notice that if x{y) is a Mobius transformation then tuq = 1 and G — 1, so the 
complete eigenstructure is unchanged but for the shift from one set of characteristic 
values to another. This is not the case for more general rational transformations, 
where other changes do happen. 

The structure of the proof of Theorem 14. II is the following. First we prove the first 
part of the Theorem (the statement on elementary divisors). This is done dividing 
the statement in three cases: 

1. xo e F and yo G F; 

2. xo G F and yo = 00; 

3. Xq = 00. 

We first prove that the statement is true for case 1, then show that this implies that 
it is true for case 2. The validity of cases 1 and 2 implies case 3. 

Finally, we prove the second part of the Theorem (the statement on minimal 
indices) with a constructive proof: we build a minimal basis of Q{y) given a minimal 
basis of P{x), and vice versa. 

5. Proof of Theorem I4.lt elementary divisors. The proof relies on the 
following Lemma, whose statement generalises [7, Proposition 11.1]. The proof of the 
Lemma and more details are given in Section [51 

Lemma 5.1. Let P{x) e (F[j;])™^p and let Q{x) = A{x)P{x)B{x) where A{x) € 
Afm(F[a;]) and B[x) G Mp(F[x]) are both regular, and suppose that G F is neither 
a root o/detA(a;) G ¥[x\ nor a root of Aet B{x) G F[a;]. Then P{x) and Q{x) have 
the same elementary divisors associated with xq. 

5.1. Case 1. Define ly := mm{m,p), and let P{x) — A{x)T{x)B{x) where A{x) 
and B(x) are unimodular polynomial matrices, T(x) —: diag((5i(a:), . . . ,d^{x)) is the 
Smith form of P{x), and Si(x) are its invariant polynomials. Let now Q{y) := 
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<i?(A(x))$(T(x))$(i?(a;)). Clearly, Q{y) and Q{y) differ only for a multiplicative 
factor [d{y)]^, A G N; moreover, both det^{A{x)) and det$(i?(x)) are nonzero when- 
ever d{y) Of- Notice that, if xq is finite, then for any j/o G Tx„ there must hold 
d{yo) ^ Of (Corollary 13.71) . Therefore, Lemma [5?T] implies that Q{y), Q{y) and S{y) 
have the same elementary divisors corresponding to j/o- 

Unfortunately, $(T(x)) may not be the Smith form of Q{y), because neither 
nor are necessarily unimodular and also because ^{Si{x)) may not 

be monic. Nevertheless, it has the form diag{[d{y)]^^ 5i{y) , . . . , [d{y)]''''S^{y)), where 
ki>k2>--->ky and 5i{y) := ^{Si{x)). From Corollarv 13.71 5i{y) and d{y) cannot 
share common roots. To reduce S{y) into a Smith form, we proceed by steps working 
on 2 X 2 principal submatrices. 



In each step, we consider the submatrix 



where 7 := fc.; and 



'idiyWS^iy) 
. [d{y)]-fS,{y) 

kj, with i < j. li "/ — (j), then do nothing; if 7 > (/), premultiply the submatrix by 
l)q{y) iF-bly)qiy)\ and postmultiply it by [^{^J [a{ylp-<p\, where q{y) = 5j{y)lk{y) 
while a{y) and h{y) are such that a{y)[d{y)]'^5i{y)+h{y)[d{y)]'f'6j{y) = [d{y)Y'5i{y)] the 
existence of two such polynomials is guaranteed by Bezout's lemma, since [d(ij)\'^ 5i{y) 
is the greatest common divisor of [diy)]^ 5i{y) and [d{y)Y'5j{y). It is easy to check 
that both matrices are unimodular, and that the result of the matrix multiplications 

[d{y)\ s,(y) ^ Hence, by subsequent applications of this algorithm and 
l"C?/)J "j \y) J 

after having defined a unimodular diagonal matrix A G F'^^'^ chosen in such a way 
that the invariant polynomials of S{y) are monic, it is possible to conclude that the 
Smith form of <^{T{x)) is either S{y) := A • S{y) or S{y) S{y) ■ A (whichever of 
the two products makes sense, depending on whether m < p or not), where 

S{y) = 61a.g{[d{y)f-5M. • ■ • , [d{y)P 5rn{y)) ■ 

Thus, the ith invariant polynomial of P{x) has a root of multiplicity t,i at if and 
only if the ith invariant polynomial of Q{y) has a root of multiplicity mp^i at y^ G Tx^ ■ 

5.2. Case 2. By definition, the infinite elementary divisors for a given poly- 
nomial matrix are the elementary divisors corresponding to zero of the reversal of 
such polynomial matrix. Therefore, in order to prove Theorem 14.11 for the case 
of t/o = 00, we have to analyse the polynomial matrix Z{y) := Rev(-gc-)(5(y) = 
y^'^[d{y~^)]sp{x{y~^)), and find out what its relation to P{x) is, with particular 
emphasis to its elementary divisors corresponding to j/o = Op. Recalling Remark 13. 5[ 
notice that there are two distinct subcases for which 00 G T^g for a finite xq G F. We 
will consider them separately. 

5.2.1. Subcase 2.1: N = D = G, xq ^ nad^K We get x{y-^) = ^Zd(l) and 



y'^d{y ^) = Kevd{y); therefore Z{y) = [Rev(j(i/)]^P( |^°^"^^| ). This means that we 



• 6^" -^yy I — Rovd(y) 

- Rcvn(y) ' 
' Rcvti{y) - 
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can prove analogous results for Z{y) just by considering this time the new rational 
transformation y ^ x — . From Remark 13.51 Op is a root of multiplicity 

G — S for the equation Revn(?/) — xoHevd^y); moreover, since we took the reversal 
with respect to the degree (or also because of Corollarv I3.7p . Of cannot be a root 
of Kevd{y). Therefore, following the proof given above, one can state that P{x) has 
{x — xoY^ , ■ . ■ ,{x — xoY^ as elementary divisors corresponding to xq if and only if 
Z{y) has the j elementary divisors . . . ^ y(G-S)ij corresponding to Of- The 

thesis follows immediately. 

5.2.2. Subcase 2.2: N < D ^ G, xq = Of- This time, we can write x{y^^) = 
'''neXf'^ ^^^'i ^(y) = [Revd(y)]^ P{y^^^). It is therefore sufficient to 
consider the transformation y ^ x ~ y^~'^ Rcvdfa) • 

In fact, notice that Of is a solution of multiplicity G — N for the equation 
y^~^Revn{y) = (Of is neither a root of Revn(y) nor a root of Revd{y), because 
Revn(OF) — ^ Of and Revd(OF) = do ^ Of). Thus, P{x) has the j elementary di- 
visors , . . . , corresponding to Of if and only if Z{y) has the j elementary divisors 
y{G-N)ii^ . . . , corresponding to Of, and the thesis follows. 



5.3. Case 3. By definition, the infinite elementary divisors of P{x) are the ele- 
mentary divisors corresponding to the characteristic value Of for R{x) :— RevgP(a;) = 
xfP(a;"i). But let *g,„(y),d(a) = g ,d(y) ,n(v) and U{y) = *g(i?(x)), that is to say 
U{y) — [niy)]^ R{^^^) . A simple calculation gives 

U{y) = [niy)]^[^]^P{[^]-') = [d{y)]^Pi^) = <f,(P(2/)) = Qiy). 

niy) n(y) d(y) 



One can therefore follow the proof as in the previous Subsections, but starting 
from R{x) and using a different transformation (notice that the equation d{y) — Op 
defines both Too for the old transformation and Tqj. for the new transformation). 

6. Proof of Theorem [471} minimal indices. We shall only prove the theorem 
for right minimal indices. The proof for left minimal indices follows from the proof 
for right minimal indices and from the fact that $ and the operation of transposition 
commute, that is ^g{P^{x)) = {^g{P{x))Y V P{x) e {¥[x])"'''p. 

6.1. =>. Let dimkerP(a;) = s, and V{x) — [vi{x), . . . ,Vs{x)] be a minimal ba- 
sis for kerP(a;), with minimal indices /3i :— degw^ Vi — l,...,s and order B := 
J2i=i Pi- Foi' each value of i let us define Wi{y) := (wi(a;)); we also define W{y) := 
[wi{y), ... , Ws{y)]. Clearly deg Wi{y) ^ G^i. Suppose in fact deg Wi{y) ^ Gfii, apply- 
ing Proposition 13.21 (in the case g ~ k ~ l3i), this would imply that there exists some 
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G F and some polynomial vector u(x) G (F[a;])^ such that Vi{x) — (x — xo)u{x). 
Hence, [vi{x), . . . ,{x — xo)~^Vi{x), . . . ,Vs{x)] would be a polynomial basis of order 
B — 1 for kerP(a;), leading to a contradiction. In order to prove that W{y) is a 
minimal basis for ker(5(a;) we must show that it is a basis and that it is minimal. 

Clearly Wi{y) lies in kerQ(?/) for all i. In fact, P{x)vi{x) = implies that 
Q{v)wi{y) — 0. So it is sufficient to show that W{y), considered as an element 
of (F(a;))P^^ has rank s. Notice that W{y) = V{x{y)) ■ dvAg{[d{y)Y\ . . . , [d{y)]f^^). 
A well-known property of the rank is that, if Ai — A2A3 and ^3 is square and 
regular, then rk(Ai) = rk(A2). Therefore rk(W^(y)) — Tk{V{x{y)), because the di- 
agonal matrix above is regular. Let V{x) be some regular s x s submatrix of V{x), 
which exists because rk(V^(a;)) — s. By hypothesis, det(F(a:)) ^ Of[x], which implies 
det{V{x{y))) ^ O^y). Hence s = rk{V{x{y))) = T]i{W{y)). Then W{y) is a basis. 

In order to prove that it is minimal, let us introduce the following lemma whose 
proof can be found in [2] . 

Lemma 6.1. Let V be a vector subspace ofWlxY, with dimV = s. Let H = 
[hi, ... , hg\ be a polynomial basis of order A for V and define ^i, i = 1, . . . , ( f ) to be 
the s X s minors (i.e. determinants of s x s submatrices) of H. Then the following 
statements are equivalent: 

• H is a minimal basis for V 

• The following conditions are both true: (a) GCD(^i, . . . , = 1f[2;] (b) 
maxi deg^,; — A. 

So let ^i(y) be the s x s minors of W{y). We shall prove that (a) their GCD is 
1f[j^] and (b) their maximal degree is GB = GJ2i=if^i- By Lemma ISTTl these two 
conditions imply that W{y) is minimal. Recall that Wi{y) = that is to 

say Wi{y) = [d{y)]^^Vi{x{y)). Any s x s submatrix of W{y) is therefore obtained from 
the corresponding s x s submatrix of V{x) by applying the substitution x = x{y) 
and then multiplying the ith column by [d{y)]^^ for i — l,...,s. Let us call (i{x) 
the s X s minors of V{x). From the properties of determinants we obtain the relation 

= {UUMy)f') C^ix{y)) = [diy)^Mxiy)) = $b(C.(x)). 

Now for each i let ji := deg Ci (x) and 5i :— deg [y) < maxj<^. {Nj — Dj) + DB 
where the maximum is taken over those values of j such that the jth coefficient of (y) 
is nonzero (Lemma l3.1|) . There are two cases, li N < D = G, 6i < GB, and applying 
ProDOsition l3.2l (with g ~ B), the inequality holds if and only if if [x — x)\(^i{x) , where 
X ^ \i N < D and x = ncd^^ if N = D; notice that there must be at least 
one value of i for which 6i = GB, otherwise (x — x) would be a common factor of 
all the Ci{x)j which is not possible because of Lemma (6. II Finally, if D < iV = G, 
di = ^iG + {B — ^i)D. Since V{x) is minimal we have maxi(7i) = B, which implies 
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that also in this case ma.Xi{Si) = GB. This proves condition (b). 

Notice moreover that S,i{y) = '^B{C,i{x)) = [d{y)\^~'^^^^-{C,i{x)), where the first 
and the second factor are coprime (because of Corollary 13. 7p . Let us prove the fol- 
lowing Lemma. 

Lemma 6.2. Let p,q,r ^¥[x\ withr monic. Then, GCDj-jx] (p, (/) —r if and only 
i/ GCDF[y](<i>dogp(p), $dcgq(g)) = K ■ $dcgr(f), where k e ¥ is such that k ■ ^dogr{r) 
is monic. 

Proof. Let a, (3 be two suitable elements of F and let us write the prime factor 
decompositions p — a-Y\{x — piY^ , q = (3 ■Y\{x — qi)^^ , r = Y\{x — ri)''' . Of course 
we have that {x — riY^ \r if and only if [x — PiY^ \p, {x — qiY*\q and pi — min(7ri, 9i). 
We get 3>dogp(p) = a ■ ]\{n{y) ~ pd{y)Y% ^dog^Cg) = P ■ ]\{n{y) ~ qtd{y)Y' and 
'i'dcgr('') = n("(2/) ~ i^id{y))''' ■ The thesis follows by invoking Corollary 13.71 □ 

Lemma \6l2\ implies condition (a). This follows from the equation GCDi(fi(y)) = 
GCDi([d(?/)]^~''') • GCDi($^. (^i(x))) = lw[y] ■ lw[y], where the first lf[y] comes from 
the fact that maxi(7i) = B, while the second ^wiy] comes by applying the previous 
Lemma to GCD(ei (y),..., Uv)) = GCD(GCD(. . . GCD(e2(2/), ^(y)) . • . )) and from 
the identity <I>o(1f[:e]) = l¥[y]- 

6.2. -4=. To complete the proof, suppose now that Q{y) = ^g{P{x)) for some 
P{x) e ¥[x] and that W{y) is a minimal basis for ker(3(?/), with minimal indices 
El < ••• < Es- The other implication that we proved in the previous subsection 
implies that G\ei V i, so define l3i — ^. Suppose that there exists a minimal basis 
V{x) — {vi{x), . . . , Vs{x)) for ker P{x); suppose moreover that an index G {1, . . . , s} 
exists such that degVi„ ^ jSi^. Applying the reverse implication, this would imply 
that there is a minimal basis W{y) = {lii (y), . . . , Ws{y)) for ker Q{y) whose ioth right 
minimal index is not equal to . This is absurd because every minimal basis has the 
same minimal indices. 

7. Extension to more releixed hypotheses. For the sake of convenience, 
we have so far assumed that the field F is algebraically closed. This unnecessary 
hypothesis can be dropped. To see it, assume that F is not algebraically closed and 
let K be the algebraic closure of F. Then {¥[x])"'''p C (K[a;])™^P, so we can use 
Theorem 14.11 to identify the Smith forms of P{x) and Q{y) — $g(P(a;)) over the 
polynomial rings M.[x] and IK[?/]. We can then join back elementary divisors in K[a;] 
and K[y] to form elementary divisors in ¥[x] and ¥[y]. Of course, in this case an 
elementary divisor is no more necessarily associated with a characteristic value in F. 
For instance, if F = Q, then the elementary divisor a;^ + 2 is not associated with any 
rational characteristic value, but if we consider the field of complex algebraic numbers 
K = Q then we can split it as {x — \/2i){x + \/2i) and associate it to the characteristic 
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values ±\/2i- Similarly, the other results (e.g., Lemma [6.2p that use the algebraic 
closure of F can be straightforwardly extended to a generic field F via an immersion 
into its algebraic closure K. 

8. Proof of Lemma mi Let P(a;) e F™^p[x]. IiU{x) := [ui{x), Us{x)] is 
a minimal basis for kerP(a;), we define ker^o P{x) :— span ({ui(a;o), . . . , Us{xo)}) C 
F^*. In general ker^^g P{x) is a subset of kerP(a;o). It is a proper subset when xq is a 
characteristic value of P{x), as is illustrated by the following example: let F = C and 

r a; 1 1 



. a; . 

Evaluating the polynomial at 0, we get kerP(O) = span({[l, 0, 0, 0]'^, [0, 0, 0, l]"^}). 
On the other hand, a minimal basis for kerP(x) is [l,~x,x , Of, so kero P{x) = 
span({[l, 0,0,0]^}). 

We need now to slightly modify a definition given in [7] in order to extend it to 
the case of singular and/or rectangular polynomial matrices. A polynomial vector 
v{x) £ {¥[x])P is called a root polynomial of order £ corresponding to xq for P{x) if 
the following conditions are met: 

1. Xq is a zero of order £ for P{x)v{x); 

2. v{xo) ^ ker^„ P{x). 

Observe that v{xo) S ker^^j, P{x) 4^ 3 w{x) e kcrF(a;) C {W{x)y : w{xo) — v{xo)- 
In fact, let w{x) = U{x)c{x) for some c{x) G (F(x))* and w{xo) = v{xo): then 
v{xo) = U{xo)c{xa) € ker^iQ P{x). Conversely, write v{xo) = U{xq)c for some c € F^ 
and notice that U{x)c e kerP(a;). Hence, condition 2. implies v{x) ^ kerP(x). 

In [71 Proposition 1.11] it is shown that given three regular polynomial matrices 
P{x), A{x), B{x) £ AIn{x), and if Xq is neither a root of detA{x) nor a root of 
deti?(x), then v{x) is a root polynomial of order £ for A(x)P{x)B{x) corresponding 
to Xq if and only if B{x)v{x) is a root polynomial of order £ corresponding to xq for 
P{x). The next Proposition generalises this result. 

Proposition 8.1. Let P{x) e F™^p[x], A{x) e i\/™(F[x]) and B{x) e Mp{¥[x]). 
Suppose that both A{xo) and B{xo), with xq S F, are full rank matrices. Then v{x) 
is a root polynomial of order £ corresponding to xq for A{x)P{x)B(x) if and only if 
B{x)v{x) is a root polynomial of order I corresponding to Xq for P{x). 

Proof. Notice that if A{xq) and B(xq) are full rank then A{x) and B{x) are 
regular. In [7], a root polynomial is defined for regular square polynomial matrices, 
so that condition 2. reduces to v{xq) ^ 0. Nevertheless, the proof given in [71 
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Proposition 1.11] for condition 1. does not actually use the regularity of P(a;), and it is 
therefore still valid when P{x) is not a regular square polynomial matrix. To complete 
the proof: v{xo) £ ker^g A{x)P{x)B{x) <^ 3wi{x) £ kei A{x)P{x)B{x) : wi{xo) = 
v{xo) <^ 3^2(2:) S kerP(a:) : W2{xo) = B{xo)v{xo) <^ B{xo)v{xo) G ker^^ P(a;). To 
build ^2(2;) from wi{x), simply put W2{x) — B{x)'Wi{x) and use the fact that A{x) is 
regular. To build wi{x) from W2{x), let {B{x))~^ be the inverse matrix (over F(a;)) 
of B{x), which exists because B(x) is regular; then, put wi{x) = {B{x))^^W2{x). □ 

Let vi{x), . . . , Vs{x) be root polynomials corresponding to xq of orders ii < ■ ■ ■ < 
is- We call them a maximal set of Xq -independent root polynomials if: 

1. they are a;o-independent, i.e. fi(a;o), . . . ,Vs{xo) are linearly independent; 

2. no (s+l)-uple of xo-independent root polynomials corresponding to a;o exists; 

3. there are no root polynomials corresponding to xq of order £ > £s] 

4. for all J = 1, . . . , s — 1, there does not exist a root polynomial Vj{x) of order 
ij > ij such that Vj (x) , Wj+i (x) , . . . ,Vs (x) are xo-independent. 

As long as det B{xo) and det A{xo) are nonzero, it is easy to check that vi{x), . . . , Vs{x) 
are a maximal set of XQ-independent root polynomials for A{x)P{x)B{x) if and only 
if B(x)vi{x), . . . , B{x)vg{x) are a maximal set of XQ-independent root polynomials for 
P{x). The next Proposition completes the proof of Lemma |5. II 

Proposition 8.2. P{x) e {¥[x])"^^p has a maximal set of xo-independent root 
polynomials, of order £1, . . . ,£s, if and only if {x ~ xqY^ , . . . , {x — XqY" are the ele- 
mentary divisors of P{x) associated with xo- 

Proof. Let S{x) be the Smith form of P{x), and recall that the inverse of a 
unimodular polynomial matrix is still a unimodular polynomial matrix [7]. Thus, in 
view of Proposition 18.11 and Theorem 12.41 it suffices to prove the thesis for S{x). If 
S{x) is the zero matrix, it has neither a root polynomial nor an elementary divisor, 
so there is nothing to prove. Otherwise, let v be the maximal value of i such that 
{S{x))ii ^ Oy[x] and for j = 1, . . . ,p let ej G (F[a;])^' be the polynomial vector such that 
{ej)i = Sij. liv < p, [e^+i, . . . , Bp] is a minimal basis for ker S{x) and, being of order 0, 
also for ker^jp S{x). Suppose that vi{x), . . . , Vs{x) is a maximal set of a; 0- independent 
root polynomials for S{x). Let fc < be the smallest index such that {vs{xo))k 7^ 0^: 
there must exist such an index because Vs{xo) ^ keixg S{x). Let {S{x))kk = ~ 
xo)''9{x), with 9{xa) ^ Of. We get {S{x)vs{x))k ^ {x - XQy9{x){v4x))k, so ^ > is. 
Actually, ^j, = £s, or would be a root polynomial of order greater than £s , which is 
absurd. Then let k' be the largest index not equal to k and such that {vs-i{xo))k' ^ Ow 
(if such and index does not exist, then Vg-iixo) is, up to a vector lying in keixo S{x), 
a multiple of and thus £s-i — £s'- in this case, replace without any loss of generality 
Vs-i{x) by a suitable linear combination of Vs-i{x) and Vs{x)). Following an argument 
similar as above, we can show that {S{x))k'k' — {x ~ xoY'^^9{x), 9{xo) 7^ Of. We 
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repeat the process until we find all the s sought elementary divisors. There cannot be 
more, otherwise dimkerS'(xo) — dimker^jo S{x) > s and it would be possible to find 
an (s + l)-uple of xg-independent root polynomials. Conversely, it is easy to check 
that Ciy-s+i, ■ ■ ■ ,61, are a maximal set of XQ-independent root polynomials. □ 

Remark 8.3. Root polynomials carry all the information on Jordan chains J^. 
Letvix) = "Y^l^liix — xoYvi, Vi € F™, be a root polynomial of order £ corresponding to 
xq for P{x). It is possible to prove that then w{y) = Yfi=o['^iy)Y^^^^[''^iy) ~ ^od{y)Yvi 
is a root polynomial of order n^oi corresponding to yo for Q{y). The latter formula 
relates the Jordan chains of Q{y) at yo to the Jordan chains of P{x) at xq. 

9. Conclusions. We have shown that if P{x) and Q{y) are polynomial matrices 
whose entries belong to the ring of univariate polynomials in x (resp. y) with coef- 
ficients in any field, and if P{x) and Q{y) are related by a rational transformation 
x{y)^ then the complete cigenstructures of Q{y) and P{x) are simply related. 
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